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We characterize the long-time projective behavior of the stochastic master equation describing a 
continuous, collective spin measurement of an atomic ensemble both analytically and numerically. 
By adding state based feedback, we show that it is possible to prepare highly entangled Dicke states 
deterministically. 
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I. INTRODUCTION 

It has long been recognized that measurement can be 
used as a non- deterministic means of preparing quan- 
tum states that are otherwise difficult to obtain. With 
projective measurements that are truly discrete in time, 
the only way an experimentalist can direct the outcome 
of the measurement is by preparing the initial state to 
make the desired result most probable. Generally, it is 
impossible to make this probability equal to one, as the 
measurement will, with some non-zero probability, result 
in other undesirable states. If the experimentalist can af- 
ford to be patient, then accepting a low efficiency is not a 
problem, but this is not always the case. In recent years, 
a theory of continuous quantum measurement has been 
developed that fundamentally changes the nature of state 
preparation via measurement Ij. When a measurement 
and the corresponding acquisition of information are suf- 
ficiently gradual, there exists a window of opportunity 
for the experimentalist to affect the outcome of the mea- 
surement by using feedback control In this paper, we 
demonstrate that it is possible to deterministically pre- 
pare highly entangled Dicke states 0,0 of an atomic spin 
ensemble by adding state based feedback to a continuous 
projective measurement. 

It has been shown that models of quantum state reduc- 
tion exist that exhibit the usual rules of projective mea- 
surement except the state reduction occurs in a continu- 
ous, stochastic manner These models are not without 
physical relevance as they are the same as those derived 
to describe the conditional evolution of atomic spin states 
under continuous quantum non-demolition (QND) mea- 
surement 0, 01 • By measuring the collective angular mo- 
mentum operator, J^, of an initially polarized coherent 
spin state via the phase shift of an off-resonant probe 
beam it ispossible to produce conditional spin squeezed 
states 1^ I3 which are of considerable interest for appli- 
cations in quantum information processing and precision 
metrology [Him. 

In these models, the reduction in variance that initially 
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leads to conditional spin squeezing is the pre-cursor of 
the projection onto a random eigenstate of at longer 
times. Figure n demonstrates the projection process for 
a single numerically simulated measurement trajectory 
[3^ . Like spin squeezed states, these Dicke states offer 
potential for quantum information applications because 
of their unique entanglement properties jl2| |. Although 
the experimental difhculties in obtaining these states via 
QND measurement or other experimental methods [l^ 
in . [T5I I are considerable, the details of the continuous 
projective process that leads to them are of fundamental 
interest. 

Whenever the measurement is sufficiently slow, an ex- 
perimentalist may steer the result by feeding back the 
measurement results in real time to a Hamiltonian pa- 
rameter. Indeed, the measurement process, as a state 
preparation process, can be made deterministic with the 
use of feedback control. Recently, we have demonstrated 
this concept by modulating a compensation magnetic 
field with the measurement record to deterministically 
prepare spin squeezed states 9f as proposed in ^f!\. This 
is just one example of the growing confluence of quantum 
measurement with classical estimation and control the- 
ory 0, 0. Other applications of quantum feedback 
include parameter estimation, metrology, and quantum 
error correction (IS. 19. 20. .21112^. 

In this paper, we focus on the long time limit of 
the QND measurement and feedback process. Just as 
spin squeezed states can be deterministically prepared at 
short times, we numerically demonstrate that individual 
Dicke states can be deterministically prepared at long 
times with the use of state based feedback While 
our proposed feedback laws are non-optimal, they demon- 
strate the adequacy of intuitive controllers with finite 
gain for directing the diffusion of the quantum state to- 
wards desirable regions of Hilbert space with unity effi- 
ciency. This is in contrast to other proposed schemes us- 
ing measurement to prepare Dicke states probabilistically 

dill. 

A more systematic approach utilizing stochas- 
tic notions of stability and convergence in the continuous 
measurement and control of a single spin is presented in 
reference 0. 

This paper is organized as follows. In Sectional we in- 
troduce the stochastic master equation which represents 
the rule for updating the system state in real time via the 
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FIG. 1: The results of a single numerical simulation [36] of the SSE, Eq. 0, with M = 1, 77 = 1, and A'' = 10 spins initially 
aligned along the a;-axis. (A) In a quantization axis perpendicular to the polarization, the level distribution of a coherent spin 
state (CSS) is Gaussian for large N. Under conditional measurement the state evolves at short times into a spin squeezed state 
and, eventually, into a random eigenstate of . (B) A map of the state's angular distribution on the Bloch sphere in spherical 
coordinates. The uncertainty in the transverse direction to the measurement axis grows until there is no information about the 
perpendicular component direction. (C) At long times, the population is at most divided among two levels that compete to be 
the final winner, which in this case appears to be m = 1. (D) All of the state information is obtained by properly filtering the 
noisy photocurrent. 



incoming measurement record. Here we discuss the var- 
ious representations of the dynamics in both the short 
and long time limits. Section IIIII describes the proba- 
bilistic preparation of Dicke states via observation alone. 
The numerical demonstration of the open loop projec- 
tion process reveals statistical features that clarify the 
details of the projection. Feedback is added to the pro- 
cedure in Section IIVI where we show that state based 
control allows one to prepare the same Dicke state deter- 
ministically on every measurement. Finally, in Section 
Ivl we discuss future directions and imminent challenges 
regarding quantum state preparation via measurement 
and control. 



II. REPRESENTATIONS OF THE 
CONDITIONAL EVOLUTION 

The physical system we will consider is an ensemble 
of N spin-1/2 particles contained within a cavity and 
interacting with a far off resonant single field mode. We 
will denote the conditional state of the spin ensemble 
as p{t) and the homodyne measurement record of the 
output as y{t). The stochastic master equation (SME) 



describing the conditional evolution is @, 01 

dp{t) = -i[H{t),p{t)]dt + V[y/MJ;:]p{t)dt 
+ V^W[VMJ.]p(<) (2^/M^[yit)dt ~ {.h)dt]^ (1) 

where H{t) = ^Jyb{t) is the control Hamiltonian that we 
will allow ourselves (without feedback h{t) = 0), 7 is the 
gyromagnetic ratio, M is the probe parameter dependent 
measurement rate, and 

V[c]p = cpc' - [c'cp + pc'c)/2 (2) 
H[c]p = cp pc^ - Tr[(c + ct)p]p (3) 

The (scaled) difference photocurrent is represented as 

y(t)dt = {J:,){t)dt + dW{t)/2y/M^ (4) 

The stochastic quantity dW{t) = 2^/ M-q{y{t)dt — 
{Jz)[t)dt) is a Wiener increment and dW{t)/dt is a Gaus- 
sian white noise that can be identified with the shot-noise 
of the homodyne local oscillator. (See [25I |2^ for an 
introduction to stochastic differential equations (SDE).) 
The sensitivity of the photodetection per VHz is repre- 
sented by l/2v^M?7, where the quantity 77 G [0, 1] repre- 
sents the quantum efficiency of the detection. If ry = 0, we 
are essentially ignoring the measurement result and the 
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conditional SME becomes a deterministic unconditional 
master equation. If = 1, the detectors are maximally 
efficient. In this latter case, the conditioned state will 
remain pure for the entire measurement, thus we can use 
a state vector description, and the SME can be replaced 
with a stochastic Schrodinger equation (SSE) 

dm)) = i-iHit) - M{J, - {J,){t))V2)\^{t))dt 

+VM{J, - {J,){t))m))dW{t) (5) 

This SSE was considered in where the motivation was 
more abstract and less concerned with the experimen- 
tal filtering perspective presented here. We emphasize 
that the SME/SSE is physically derived and is an ex- 
plicit function of a measured photocurrent variable y{t), 
through which the randomness enters. The states are 
considered as states of knowledge and, in practice, an ex- 
perimentalist updates the description of the system, p(t) 
(Figure njV-C) , as the measurement results, y{t) (Figure 
^3), arrive in time. 

Even though the light field is dispersively coupled to 
each atom, there will still be some absorption of light and 
hence decoherence for which the above master equation 
does not account. If the cavity coupling is strong enough, 
and extraneous noise sources are sufficiently reduced, the 
measurement rate will be faster than this decoherence 
and we can consider the long time limit of the master 
equation 0, Q- For free space measurements H^, H^, 
e.g. free space Faraday rotation or homodyne detection 
of one polarization component, the above SME will be 
approximately valid in the short time limit only. 

A. Hilbert space, coherent spin states, and Dicke 
states 

Under certain idealizations, we can considerably re- 
duce the size of the Hilbert space needed to describe the 
conditionally measured ensemble. Throughout this pa- 
per, the initial state p{0) will be made equal to a coher- 
ent spin state (CSS) polarized along an arbitrary direc- 
tion For example, a CSS pointing along the z-axis is 
denoted | tiT2 • • • '\n)z and all others can be prepared 
by rotating this state with the angular momentum oper- 
ators Ji, with i € {x, y, z}. A CSS, typically obtained via 
a dissipative optical pumping process, is an eigenstate of 
with maximal eigenvalue J( J -I- 1), where J — N/2. 
Because the SME works under the QND approximation 
of negligible absorption (i.e. the large detuning dispersive 
limit), no angular momentum will be exchanged between 
the probe beam and the ensemble. The measurement and 
possible field rotations represent the only allowed dynam- 
ics, thus the state will maintain maximal (J^) over the 
course of the measurement. 

The Dicke states are defined as the states |/, m) that 
are simultaneous eigenstates of both and Jz 

Jz\l,m) = m\l,m) (6) 
j2|Z,m) = l{l + l)\l,m) (7) 



where 

\m\<l< J = N/2 (8) 

Under the above approximations, we can neglect any 
state with / ^ J. We then shorten the labelling of our 
complete basis from | J, m) to | m) so that 

Jz\m) ~ m\m) (9) 
J2|m) = J{J+l)\m) (10) 

where m S {-N/2, -N/2 + 1, N/2 - 1, N/2}. 

When the physical evolution is such that the |m) states 
remain complete, we can limit ourselves to a density ma- 
trix of size {N -I- 1) X {N + 1) rather than the full size 
2^ X 2^. This reduced space is referred to as the sym- 
metric subspace, as its states are invariant to particle 
exchange |22, |30|. For the case of two spins, the sym- 
metric subspace contains the triplet states, but not the 
singlet. States contained within the symmetric subspace 
can be described as a pseudo-spin of size J = N/2. 

In the z basis, the extremal values of m, ±N/2, are 
simply the coherent spin states pointing along the z-axis 

\m = +N/2) = |tiT2---Tiv) (11) 
\m = -N/2) = Uii2---M (12) 

In terms of the constituent spins, these states are obvi- 
ously unentangled. In contrast, consider the state with 
m = (for iV even) 

|m = 0) - CI],P,(| Ti • • • U/2iN/2+i ■ ■ ■ In)) (13) 

where the represent all permutations of the spins and 
C is a normalization constant. This state is highly en- 
tangled in a way that is robust to particle loss 12]. Even 
though the expectation values (Ji) vanish for this state, 
it still has maximal eigenvalue. Loosely, this state 
represents a state of knowledge where the length of the 
spin vector is known and the z component is known to 
be zero, but the direction of the spin vector in the x — y 
plane is completely indeterminate. Similarly, the entan- 
gled states with < |m| < N/2 can be imagined as living 
on cones aligned along the z-axis with projection m. The 
loss of pointing angle information from the measurement 
process is diagrammed in Figure ^3. 

Along with their unique entanglement and uncertainty 
properties, Dicke states are also of interest for the impor- 
tant role they play in descriptions of collective radiation 
processes and for their potential role in quantum in- 
formation processing tasks ,13i . ^14 iSL] . 

B. Short time limit 

Even when working within the symmetric subspace, 
for a large number of spins the size of p{t) may be too 
unwieldy for computational efficiency. Because it is of- 
ten desirable to update our state description in real-time 
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(e.g., for optimal feedback procedures), finding simple 
but sufficient descriptors is of considerable importance. 

We can derive a reduced model by employing a moment 
expansion for the observable of interest. Extracting the 
conditional expectation values of the first two moments 
of Jz from the SME gives the following scalar stochastic 
differential equations: 



d{J.){t) 
d{AJ!){t) 



j{J,)it)bit)dt 

+2v/M^(Ajf)(t) dW{t) 
-AM7^{^Jlf{t) dt 
-ij{[AJlJy]){t)b{t)dt 



(14) 



^{AJ^){t)dW{t) (15) 



Note that these equations are not closed because higher 
order moments couple to them. 

At short times, t <C l/ijM, we can make this set of 
equations closed with the following approximations. If 
the spins are initially fully polarized along x then, by 
using the evolution equation for the x component, we 
can show {Jx){t) ~ Jexp[— Making the Gaus- 
sian approximation at short times, the third order terms 
(AJ3) and -ij{[AJ^, Jy]){t)b{t) can be neglected. The 
Holstein-Primakoff transformation makes it possible to 
derive this Gaussian approximation as an expansion in 
1/J ■ Both of the removed terms can be shown to 
be approximately l/J\/j smaller than the retained non- 
linear term. Thus we can approximate the optimal solu- 
tion with 



d{Jz)sit) = -fJexp[~Mt/2]b{t)dt 

-f2v/M^(Aj2),(t) dWsit) 
it) = -AMr^{AJl)l{t)dt 



d{AJ: 



(16) 
(17) 



where the s subscript denotes the short time solution 
and dWs{t) = 2^/Wq[y{t)dt - {Jz)s{t)dt\. Also b{t) is 
assumed to be of a form that keeps the total state nearly 
pointing along x. The differential equation for the vari- 
ance {AJl)s{t) is now deterministic. It can be solved to 
give 



(AJf)s(i) = 



(AJ,2)(0) 



1 + 4M?7(AJf)(0)t 



(18) 



The deterministically shrinking value of {AJ'^)s{t) repre- 
sents the squeezing about the initially fluctuating value 
of {Jz)s{t) as shown in the first two frames of Figurenj^- 
B. If feedback is added, then the value of {Jz)s{'t) can 
be zeroed via Larmor precession due to a control field 
along y and the same centered spin squeezed state can 
be prepared on every trial d, ■ 

The resulting spin squeezed states can be used in sub- 
sequent precision measurements 10, llj. It is also worth 
pointing out that a precision measurement can be per- 
formed during the production of the conditional spin 
squeezing. For example, we have shown that the by prop- 
erly estimating both the spin state and an unknown clas- 
sical field simultaneously with continuous measurement 



and Kalman filtering techniques, the field estimation can 
be improved over conventional limits by the presence of 
the simultaneous squeezing 0, 0| . 



C. Long time limit 

The approximations made in the previous section are 
no longer valid at times t 3> l/r/M. The third-order 
terms become non-negligible at long times, hence the 
variance becomes stochastic. Subsequently, other high 
order moments couple into the problem and we are forced 
to consider the stochastic differential equation for each. 
Eventually, any finite numbered moment description is 
no longer useful and it initially appears that we must re- 
sort back to the full symmetric density matrix and the 
SME, Eq. as our primary description. 

Fortunately, we can take another approach and de- 
scribe the state in terms of other sufficient statistics. 
Without a field, the only statistic of the photocurrcnt 
needed to describe the state at time t is its integral, 
J*y{s)ds (see Appendix IXI or [s^l. Knowing that the 
state is only a function of this variable and the initial 
state (prior information) makes the experimental design 
of a real time estimator experimentally convenient. For 
example, we could use an analog integrator to create this 
sufficient statistic from the raw photocurrcnt, then feed 
it into a possibly non-linear device (like an FPGA :2S|) 
to perform the estimation. 

With the integrated photocurrcnt and the initial state 



(19) 



■m— — J 



we can calculate (see Appendix the conditional ex- 
pectation value of any power of Jz with the expression 



+AmMr] I y{s) ds 



(20) 



where p{t) is the unnormalized density matrix, and set- 
ting k — Q represents its trace, so 



J,*)(0=Tr[J,^p(t)]/Tr[J,V(i)] 



(21) 



Consider the case when the system starts in the x- 
polarized spin-coherent state. To very good approxima- 
tion (with reasonably large J) we can write for this state 
in the 2;-basis 



(X exp 



m 



(22) 



Using these coefficients, we now have the rule for map- 
ping the photocurrcnt to the expectation of Jz 



(J.)(t)=Tr[Jip(t)]/Tr[Jl'p(t)] 



(23) 
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Other than the minor approximation of the initial co- 
efficients, using this estimate is essentially the same as 
using solution to the full SME, so we do not give it a 
new subscript. 

To simplify further, we can change the sums to inte- 
grals, giving 



Tr[m'=p(t)] 



dm (24) 



with 



1 /■* 

A = - + 2Mr]t B = 2Mtj / y{s) ds (25) 

J Jo 



This approximation produces an estimate 



-Am^+2Bm 



dm 



(26) 



that performs sub-optimally when the distribution of 
states becomes very narrow at long times. Interestingly, 
the integral approximation here numerically appears to 
give the same estimate as the one derived previously for 
short times when no field is present, i.e. 



(JMt) = {Jz)s{t) 



(27) 



This is not entirely surprising as both of these estimators 
ignore the discreteness of the Dicke levels. Also, at long 
times, it turns out that both of these estimates appear to 
be numerically equivalent to the simplest of all estimates: 
averaging the photocurrent. In other words, one simple 
and intuitive approximation to the optimal {Jz){t) would 
be 



Ioy{s)ds 



(28) 



which is an estimate one might guess from the form of the 
photocurrent, Eq. From simulation, it appears that 
this estimate is the same as both {Jz)i{t) and {Jz)s{t) 
for t ^ l/rjM. Despite the non-optimality of these sim- 
ple estimators, they perform well enough to resolve the 
discretization of the Dicke levels at long times. 

Unfortunately, the addition of a feedback field makes 
these simplified estimators inadequate at long times, and 
deriving simple reduced models with a field present is 
difficult, thus forcing us to use the full SME in our state 
based controller. Despite this difficulty, during our sub- 
sequent feedback analysis we assume sufficient control 
bandwidth that the SME can be evolved by the observer 
in real time. 



III. MEASUREMENT EVOLUTION WITHOUT 
FEEDBACK 

In this section, our goal is to describe how the esti- 
mates of the last section probabilistically evolve at long 



times into Dicke states via observation alone. First, we 
discuss steady state and statistical properties of the SME, 
Eq. Then, we examine the unconditional dynamical 
solution with 77 which gives the average state prepa- 
ration behavior when 77 7^ 0. We then consider in detail 
how individual trajectories behave when i] ^ 0. Finally, 
we discuss the performance of the non-optimal estimators 
relative to the optimal projective estimator. 



A. Steady states of the SME and martingale 
properties 

The fact that the SME eventually prepares eigenstates 
of Jz is rather intuitive from a projection postulate per- 
spective because Jz is the quantity being measured. If 
we insert the pure Dicke state, p = |m)(77i|, into the SME 
with no Hamiltonian (or only a field along z) , we find that 
it is a steady state, dp = 0, no matter what happens with 
the subsequent measurement record. Of course, this does 
not yet prove that the state will eventually be obtained, 
as we have not discussed the stability of attractors in 
stochastic systems. 

Without a field present, the SME has several conve- 
nient properties. First of all, from the evolution equation 
for the variance notice that the variance is a stochastic 
process that decreases on average. In fact it is a super- 
martingale, in that for times s < t we have 



E,[(Ajf)(t)] < (AJ|)(.) 



(29) 



where the notation E[a;(t)] denotes the average of the 
stochastic variable x(t) at time t and the s subscript 
represents conditional expectation given a particular 
stochastic trajectory up to the time s. Additionally, it 
can be shown 5] that the average variance obeys the 
equation 



E[(Aj2)(t)] = 



(AJ.^)(0) 



I-f 4Mry(AJ2)(0)(t + C(<)) 



(30) 



where 



m 



* E[{{AJ^){s)-E[{AJ^)is)]y 
E[(AJ2)(.)]2 



■ ds 



> 



(31) 



A more explicit solution of ^{t) is not necessarily needed 
as its positivity ensures that (AJ^)(i) stochastically ap- 
proaches zero. This implies that a Dicke state is eventu- 
ally prepared. The numerical simulation of Figures^and 
13 demonstrates this behavior for an initially x polarized 
state. As expected, E[(AJ^)(t)] in Figure[5J^ appears to 
be less than the short time solution (AJ^)s(t), Eq. I|18(l . 
at long times. 

Other useful properties of the stochastic evolution are 
evident from the moment equations. For example, we 
can show that 



d{j:) 



{j:){Jz))dW{t) (32) 



012345 012345 -5 5 

time time ( J^^ ) (t) 



FIG. 2: Many open-loop moment trajectories [36] of the SSE, Eq. The trajectory of Figure Q is darkened. (A) At short 
times, the evolution of the variance (shown on a log scale) is deterministic and given by (Ajf )s(t). At long times, the variances 
become stochastic but bounded (above by 1/4 and below by exp[—2{Mt — l)]/4). The average of all 10, 000 trajectories (only 10 
are shown) gives E[(Ajf )(t)]. (B) The projective nature of the measurement is made clear by the evolution of 100 trajectories 
of {Jz){t)- The distribution of the final results is given by the first histogram of Figure UK- (C) The evolution of the 100 
trajectories all starting in an a;-polarized CSS. When rj = 1, certain regions of Hilbert space are forbidden by the evolution. 



for integer n, hence 

dm:)] = (33) 

and for times s < t we have the martingale condition 

E,[(jr)(i)] = (jr)(s) (34) 

This equation for n — 1 gives us the useful identity 

E[(J,)(t)(J,)(.)]=E[(J.)(.)^] (35) 

for s < t. Also, we can re-writc the expression for n = 2 
as 

E,[(J,)(t)2 + (Ajf)(t)] = (J,)(,s)2 + (AJ2)(.) (36) 

This implies a sort of conservation of uncertainty as the 
diffusion in the mean, shown in Figure ^3, makes up for 
the decreasing value of the variance. 

B. ri = 

It is insightful to examine the behavior of the master 
equation with 77 = which corresponds to ignoring the 
measurement results and turns the SME Eq. ^ into a de- 
terministic unconditional master equation. We continue 
to consider only those initial states that are polarized. 
This is because these states are experimentally accessi- 
ble (via optical pumping) and provide some degree of 
selectivity for the final prepared state. To see this, let us 
consider a spin-1/2 ensemble polarized in the x — z plane, 
making angle 9 with the positive z-axis, such that 

(J,)(0) = sm{9)N/2 (37) 

(J,)(0) = 

(J,)(0) = cos{9)N/2 

(AJJ)(0) = cos' {9)N/ A 

(AJ2)(0) = iV/4 

(AJ2)(0) = sin\9)N/4 



Solving the unconditional moment equations, and la- 
belling them with u subscripts, we get 

(J^)„(t) = sin(6l)exp(-Aft/2)iV/2 (38) 
{Jy)u{t) = 
(J,)„(t) = cos{9)N/2 
(AjJ)„(i) = si-a' i9)[N' ~ N ~2N'exp{-Mt) 
+ {N' - N) exp(-2Mt)]/8 + N/i 
sm''{9){N'' - N)/8 + N/4: 
(Aj2)„(t) = sm'{9)[N'~N 

+ {N - N') exp(-2Aft)]/8 + N/i 
sm'{9){N' - N)/8 + N/4: 
(Aj2)„(i) = sin2(0)7V/4 

Note that, because the unconditional solutions represent 
the average of the conditional solution, i.e. Pu{t) — 
F,[p{t)], we have 

E[(J,)(i)] = {J,)u{t) = (J.)(0) = cos{9)N/2 (39) 

This also follows from the martingale condition for 
{Jz){t). From the martingale condition for {Jl){t) we 
get 

E[((J,)(t) - E[(J,)(t)])2] 

= (Aj|)(0)-E[(Aj|)(t)] (40) 
(AJ2)(0) =sn?{9)N/4 

Thus, when < 77 < 1, we expect the final random condi- 
tional Dicke state on a given trial to fall within the initial 
z distribution. Given 9, the distribution will have spread 
I sin(6')|ViV/2 about the value cos(6')A^/2. Although the 
final state is generally random, starting with a polarized 
state clearly gives us some degree of selectivity for the 
final Dicke state because ^/N <^ N. 
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C. < ?7 < 1 

When 77 ^ 0, the measurement record is used to con- 
dition the state, and we can determine which Dicke state 
the system diffuses into. Given the task of preparing 
the state |md), the above analysis suggests the follow- 
ing experimental procedure. First, polarize the ensemble 
(via optical pumping) into an unentangled coherent state 
along any direction. Then rotate the spin vector (with a 
magnetic field) so that the z component is approximately 
equal to Tfid- Finally, continuously measure z until a time 
t ^ l/rjM. The final estimate will be a random Dicke 
state in the neighborhood of m^. When the trial is re- 
peated, the final states will make up a distribution de- 
scribed by the initial moments of Jz ((Jz)(0), (AJ^)(0), 
...). To reduce the effects of stray field fluctuations and 
gradients, a strong holding field could be applied along 
the z-axis. Because this Hamiltonian commutes with the 
observable Jz, the final open-loop measurement results 
would be unchanged. 

This process (with zero field) is shown schematically 
in Figure n for md — where the initial state is polarized 
along X. Because (Jz)(0) — 0, the final state with the 
highest probability is the entangled Dicke state md = 
0. In contrast, if (Jz)(0) = J the state would start in 
an unentangled CSS polarized along z and would not 
subsequently evolve. 

One way of characterizing how close the state is to a 
Dicke state is through the variance, (AJ^)(t). Figure[2K 
displays many trajectories for the variance as a function 
of time. For times t <^ 1/rjM the variance is approxi- 
mately deterministic and obeys the short time solution 
of Eq. l(TH|l . During this period, the mean {Jx){t) is de- 
creasing at rate M/2. Before this mean has completely 
disappeared, a conditional spin squeezed state is created. 
However, for larger times the mean and the variance 
stochastically approach zero, and the state, while still 
entangled, no longer satisfies the spin-squeezing criterion 

There are several features to notice about the approach 
to a Dicke state that are evident in Figures ^ and [3 The 
variance at time t = 1/rjM is already of order unity. 
Thus, at this point, only a few neighboring m levels con- 
tain any population, as can be seen in Figure Also, 
it can be numerically shown that, for x-polarized ini- 
tial states, the diffusion of the variance at long times 
i ^ 1/ rjM is bounded above and below by 

exph2(77A«- l)]/4 < (Aj2)(t) < 1/4 (41) 

which is evident from Figure These facts indicate 
that the population is divided among at most two levels 
at long times which 'compete' to be the final winner. If 
we assume that only two neighboring levels are occupied 
and apply the SSE (with 77 = 1), the probability, p, to be 
in one level obeys the stochastic equation 



and the variance takes the form {AJ'^){t) = p{l — p). 
As simple as it looks, this stochastic differential equation 
(SDE) is not analytically solvable jlEll^- The maximum 
variance is 1/4 and it can be shown that, for p = 1 — 
e, with e small, the lower bound is of the exponential 
form stated above, so the two-level assumption seems to 
be a good one. The fact that occupied Hilbert space 
becomes small at long times is also evident in Figure 
where the allowed states are seen to be excluded from 
certain regions when rj = 1. The arc- like boundaries of 
the forbidden space are where the two level competition 
occurs. 

In practice, an experimentalist does not always have 
an infinite amount of time to prepare a state. Eventu- 
ally spontaneous emission and other decoherence effects 
will destroy the dispersive QND approximation that the 
present analysis is based upon. Suppose our task was to 
prepare a Dicke state with, on average, a desired uncer- 
tainty, {AJ^)d <C 1, such that one level was distinguish- 
able from the next. From Eq. l|5n|) . we see that the time 
that it would take to do this on average is given by 



td 



1 



1 



(AJ|), (AJ|)(0) 



/4Mr) 



(43) 



Thus if {AJ^)d <C 1 is our goal, then td is how long 
the state must remain coherent. The larger (AJ^)(0) 
is the more entangled the final states are likely to be 
(m « 0) 0, hence, by Eq. the longer it takes to 

prepare the state for a given {AJ^)d- Hence, we arrive 
at the intuitively satisfying conclusion that conditional 
measurement produces entangled states more slowly than 
unentangled states. Of course, Eq. (|43l) is an average 
performance limit. In a best case scenario, the variance 
would attain the lower bound of Eq. (|41|1 where the state 
reduction happens exponentially fast. 



D. Performance of sub-optimal estimators 

Now we consider the performance of the sub-optimal 
estimators discussed previously, in particular the current 
average {Jz)a{t) of Eq. (|28|l . It makes sense to associate 
the overall 'error' of this estimator, denoted Va, to be 
the average squared distance of the estimator from the 
optimal estimator plus the average uncertainty of the op- 
timal estimator itself, E[(AJ^)(t)]. Using the martingale 
properties of the optimal estimate and the definition of 
the photocurrent gives this quantity as 



Va 



E[((J.) 
1 

AM-qt 



.„(t)-(J.)(t))2]+E[(Ajf)(t)] 



(44) 



dp = -2Mp{l-p) dW{t) 



(42) 



This is just the error in estimating a constant masked 
by additive white noise with the same signal to noise 
ratio 01 . The optimal estimator is better than this sub- 
optimal estimator at long times only through the quan- 
tity ^{t), Eq. 
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In the open-loop experimental procedure described at 
the beginning of the last section, the above observation 
indicates that we can replace the optimal estimator with 
the photocurrent average and still resolve the projec- 
tive behavior (given sufficient elimination of extraneous 
noise). The price paid for the simplicity of the averag- 
ing estimator is that it converges more slowly and it only 
works when a field is not present (hence without control) . 



IV. CLOSED LOOP EVOLUTION 

The primary problem with the open-loop state prepa- 
ration scheme (and other approaches 13, 14, 15j) is that 
it is probabilistic. For a single measurement, there exists 
some degree of control, by adjusting the initial angle of 
rotation 0, but the final state is a priori unpredictable 
within the variance of the initial state. In this section, 
we show that the state preparation can be made deter- 
ministic with the use of feedback. Just as the control 
scheme of 0, 0] produces deterministically centered spin 
squeezed states, we present a simple feedback controller 
that will prepare the same desired Dicke state (particu- 
larly md = 0) on every measurement trial. 

We choose to work with y-axis magnetic field actua- 
tor corresponding to the Hamiltonian, H{t) = ^b(t)Jy. 
If the CSS initial state begins in the x — z plane this 
will ensure that the vector {J){t) remains in this plane. 
This actuator is natural for the control of spin-squeezed 
states at short times, where the linear moments of (J) (t) 
are large and allow intuitive rotation of the spin vector. 
However, at long times the field will mostly be affecting 
non-linear terms in the moment expansion and the dy- 
namics are less intuitive as can be seen by the structure 
near the z-axis in Figure |2l3. Still, we continue to give 
ourselves only these rotations to work with as they are 
the most experimentally accessible actuation variable. 

In principle, the fact that Dicke states can be prepared 
deterministically with feedback should not be surpris- 
ing. Given the aforementioned characteristics of the non- 
controlled measurement one could imagine preparing a 
particular state by alternating measurement and control 
periods. For example, an initial measurement (lasting 
for a time At <^ 1 /r]M) would determine the fluctuation 
of (Jz) while the uncertainty (AJ^) simultaneously de- 
creased (on average). Then the measurement would be 
turned off and the state would be rotated with a con- 
trol field to 'zero' the conditional quantity (Jz) — rud (if 
preparing Inid) ) • The process of alternating measurement 
and control could then be repeated and would eventually 
clamp down on the desired state. Notice that, unlike the 
preparation of spin squeezed states |^ , this procedure 
could not be performed with a single measurement and 
control cycle. In other words, if we measure for a time 
t ^ l/rjM, and prepare a probabilistic Dicke state, then 
a single post-measurement rotation cannot prepare a dif- 
ferent desired Dicke state in the same basis. 

With this intuitive picture in mind, now consider the 



continuous limit of this process, where the measurement 
and control are performed simultaneously. We wish to 
find a mapping from the photocurrent history to the con- 
trol field that prepares our state of interest in a satis- 
factory manner on every trial. For simplicity, we work 
with rj — 1 and use the SSE of Eq. for all simu- 
lations [36]. In selecting a controller, we could choose 
one of several strategies, including either direct current 
feedback or a feedback rule based on the state (i.e., what 
has been called Markovian and Bayesian feedback respec- 
tively While direct current feedback possesses 
certain advantages, mainly simplicity that allows practi- 
cal implementation, and is capable of working adequately 
at short times, any constant gain procedure would never 
prepare a Dicke state with confidence. If the current is 
directly fed back, a finite amount of noise will unneces- 
sarily drive the system away from its target, even if the 
state starts there. Of course the gain could be ramped 
to zero in time, but unlike the short time case, it is not 
clear how to tailor the gain intelligently. 

Another alternative would be to prepare a spin 
squeezed state with this approach and then turn off the 
feedback at some intermediate time. This would cer- 
tainly enhance the probability of obtaining a certain 
Dicke state, but the process would remain probabilistic to 
some degree. For these reasons, we continue considering 
only state based feedback, despite the fact that updating 
the state in real time is experimentally challenging. 



A. Defining a cost 

A useful first step in the design of any controller is 
to define the quantity that the ideal controller should 
minimize: the cost function. For example, consider a 
state preparation application where the controller aims 
to produce the desired target state IV'd)- In this case, one 
possible cost function is the quantity 

Uf = l- (MpIAi) > (45) 

evaluated at the stopping time, which is zero iff the fi- 
delity of the state with respect to the target is unity. 
In the current application, where we desire a final Dicke 
state \md) we wish to minimize a different quantity 

U = ((Jz) - mdf + (Ajf ) 
= T.m{m\p\m)^ (m - nidf 
> (46) 

which is zero iff p = \md){md\- Notice that U gives a 
higher penalty than Uj to states that are largely sup- 
ported by Dicke states far removed from the target. In 
general, U will evolve stochastically and we may be more 
interested in the mean behavior, denoted E[J7]. In the 
uncontrolled case, it can be shown that this quantity re- 
mains constant, E[?7(t)] = U{0). For the controlled case, 
we wish for E[J7] — > as time advances, which, because 
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012345 012345 -5 5 



time time ( J^^ ) (t) 

FIG. 3: One hundred closed-loop moment trajectories [36] of the SSE with feedback law b{t) = X{JxJz + JzJx){t)/'2 and A = 10 
chosen from numerical considerations. (A-B) If the control is successful the quantity {J^}{t) should go to zero on every trial. 
For this controller the number of successful trajectories is increased significantly (from 25 to 90 percent), but the remaining 
fraction is attracted to neighboring fixed points, causing the mean E[(jf)(f)] to saturate at a non-zero value. Although the 
successful fraction converges exponentially, the fastest converging trajectories are slower than in the open loop case. This is 
evident in (C) as the converging trajectories have visibly not yet reached {Jx) = at time t — 5. 



U > 0, implies that every trajectory approaches the tar- 
get state \md)- 

In general, the cost function could also include an in- 
tegral of the quantity U{t) instead of just the final value. 
As in classical control theory it is also practical to 
include a function of b{t) in the cost as a way of ex- 
pressing our experimental feedback gain and bandwidth 
constraints. Analytically proceeding in this way by op- 
timizing the average cost is too difficult for the current 
problem, but, with this perspective in mind, we proceed 
by proposing controllers according to related considera- 
tions. 



B. Control law 1 

Now consider the average evolution of the above cost 
function, which is given by 



dE[U{t)] = 
-27E 



2 



Because we want this function to continuously decrease, 
the right hand side should be negative at all times. If 
we have full access to the density matrix, and minimal 
feedback delay, we could use the controller 

biit) ^ A (V£^±^M{^ _ (48) 

where A is a constant positive gain factor. This law guar- 
antees that c?E[[/(i)] < 0. Still, this does not yet prove 
that t7 = is obtained because (iE[[/(t)] = for states 
other than the target state. Furthermore, even with this 
control law applied, all Dicke states remain fixed points. 

Regardless of these issues, we proceed by analyzing the 
performance of this control law numerically with — 0. 
In principle, the gain could be chosen arbitrarily large. 



Here we choose to work with a gain that is large enough 
to be effective but small enough to keep the numerical 
simulation results valid [36] . The choice of a limited gain 
is a necessity in both simulation and experiment, thus 
we wish to find a control law that works within this con- 
straint. For the parameters used in our simulation, we 
use a gain of A = 10 which produces the results shown in 
Figure 01 

In Figure we now plot the figure of merit for 
rud = 0, U{t) — {Jl){t). In open loop configuration, only 
25% of all trajectories are attracted to m = 0, whereas 
with this controller the percentage reaches 90%. Further- 
more, most of these trajectories approach the state at an 
exponential rate close to M, as indicated by the curve 
under which 80% of the trajectories lie. Interestingly, 
this is at the expense of those trajectories that in open- 
loop approached the target state at an exponential rate 
of 2M . There is a trade-off by which the control slightly 
compromises the convergence of the best case trajecto- 
ries. 

Unfortunately, because all other Dicke states are still 
fixed points of the controlled SSE and the gain is finite, a 
small fraction (10%) of trajectories are attracted to those 
states neighboring the target state. Thus this controller 
does not appear to deterministically prepare all trajecto- 
ries into the target state and the mean E[(J^)(i)] flattens 
at a level determined by the unsuccessful fraction of tri- 
als. 



C. Control law 2 

The obvious solution to the above problem is to try a 
controller that ensures the target state is the only fixed 
point of the SME/SSE. In this section we propose the 
control law 



62(0 =A((J,)(i)-md) 



(49) 
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time time ( J^^ ) (t) 

FIG. 4: One hundred closed-loop moment trajectories [36] of the SSE with feedback law b{t) — X{Jz){t) and A = 10 chosen 
from numerical considerations. (A) The average over 10,000 trajectories suggests that with this control law the mean E[(J^)(t)] 
descends to zero exponentially and the target state is deterministically prepared. (B) Despite a number of early excursions, all 
100 trajectories shown converge to the desired value of m = 0. (C) Those trajectories that do not descend to the goal directly 
(about 10 of 100) are recycled and rotated back into the attractive region of the target state. Again, the control slightly 
compromises the best-case convergence rate and the trajectories have a non-zero (but still decreasing) (Jx) a,t t — 5. 



for which the state \md) is the only fixed point. However, 
unhke bi{t) this controller lacks the x symmetry that 
ensures c?E[J7(i)] < 0. Also, while the symmetry of bi{t) 
will allow it to lock to both sides of the Bloch sphere, 
62 (i) will only lock to one side of the sphere. 

Again, we proceed by numerically analyzing the per- 
formance of this controller for rrid = 0, with the results 
displayed in Figure 01 The gain is chosen in the same 
manner as before, which leads to the same reasonable 
choice of A = 10. In Figure 0p the fundamental nature 
of the dynamics can be seen. Close to 90% of the trajecto- 
ries are directly transported towards the target state, but 
the remaining fraction 'misses' on the first pass. Instead 
of being attracted towards other fixed points though, this 
unsuccessful fraction is recycled and rotated back onto 
the positive a;-axis where they can re-attempt conver- 
gence onto the target state. These large excursions can 
be seen in Figure as well, but they do not appear 

to dominate the net fiow. The average of 10, 000 trajec- 
tories gives a quantity E[{J^){t)] which appears to expo- 
nentially descend towards zero, implying that the state 
preparation has been made deterministic. As with the 
control of bi{t) there is again a trade-off: the trajectories 
that previously descended at the exponential rate of 2M 
converge more slowly, but still exponentially. 



V. CONCLUSION 

The purpose of this paper is to demonstrate the fact 
that the process of continuous projective measurement 
can be made deterministic with a theoretically simple 
and intuitive state-based control law. In the context of 
an atomic spin ensemble, the resulting Dicke states are 
highly entangled and otherwise difficult to reliably pro- 
duce from an initially unentangled state. 

However, there is much work to be done in the general 



field of quantum state estimation and control, of which 
this is one example. In this pursuit, it is helpful to uti- 
lize and adapt methods from the developed fields of clas- 
sical stochastic estimation and control theory. In p^ . 
for example, the problem of this paper is considered for 
a single spin with greater emphasis on technical notions 
of stochastic stability and convergence. Ultimately, we 
would like to discover constructive methods for deriving 
optimal control laws given a cost function and realistic 
actuators. 

Even with an optimal control law in hand, there is no 
guarantee that experimental implementation will be pos- 
sible. Any analysis should incorporate, among other con- 
straints, non-unity detection efficiencies and finite con- 
troller resources (bandwidth, memory, etc). For exper- 
imental application of quantum feedback, the controller 
complexity needs to be reduced to the point where the de- 
lay is minimal compared to other dynamical time-scales 
[23. As in classical control, effective model reduction 
techniques are indispensable when it comes to implemen- 
tation. 

Despite these difficulties, the increasing number of 
physical systems that can be measured reliably at the 
quantum limit will surely hasten the effort to solve many 
of these technical challenges. By respecting the physical 
basis of measurement dynamics, experimentalists will be 
able to more efficiently use measurement itself, in tan- 
dem with more traditional techniques, to actuate quan- 
tum systems into desirable states. 
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APPENDIX A: SOLUTION OF THE SME 
WITHOUT A FIELD 

An explicit solution to the SME, Eq. can easily be 
found in the case H{t) — 0. First, the SME is rewritten 
as 

dp{t) = V[\/MJ^]p{t)dt + 2Mri {Jzp{t) + p{t)Jz) y{t)dt 

(Al) 

This equation, known as the unnormalized or linear SME, 
is equivalent to Eq. with the identification 

p{t) = p(t)/Tr[p(t)] (A2) 

Introducing the notation 

Qip = JzpJz 

g2P = Jlp+~pJl (A3) 
QaP = JzP + pJz 

Eq. can be written in the more suggestive form 

dp{t) = M{g^ - \g2)p{t)dt + 2Mrj Gspit) y{t)dt (A4) 

Now note that Eq. ljA4p is a linear Ito stochastic dif- 
ferential equation (SDE) for and moreover 
Si,2,3 all commute with each other in the sense that 
QiQjP = QjQiP- Such SDEs have a simple explicit so- 
lution ^ 

p{t) = cxp[(M{l^r])gi-M(l + 7^)02 /2)t 

+2Mr,g3 I y{s)ds\p{0) (A5) 
Jo 



as is easily verified by taking the time derivative of this 
expression, where care must be taken to use Ito's rule for 
the stochastic term. 

Now consider an initial pure state of the form 

J 

m— — J 

The associated initial density matrix is then 

J 

m = \m){m\^ E c„,c,vim)(m'i (a?) 

m,m' — — J 

Substituting into Eq. l|A5p gives 
J 

p{t) = E CmC*„, exp[(Af (1 - vi)mm' 

m,m' — — J 

- iM(l + ?7)(m2 + {m!f))t 

+ 2Mr]{m + m') [ y{s)ds]\m){m'\ 
Jo 

(A8) 

Hence 

J 

Tr[J^p{t)] = E "^''|cmpexp[-2A^r/m2^ 

m— — J 

+AniMr] f y{s) ds] (A9) 
Jo 

which is the result used in the text, Eq. (|2U|I . 
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